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We study the infrared behaviour of lattice 51/(3) Yang-Mills theory in Coulomb gauge in terms 
of the ghost propagator, the Coulomb potential and the transversal and the time-time component 
of the equal-time gluon propagator. In particular, we focus on the Gribov problem and its impact 
on the observables. We observe that the simulated annealing method is advantageous for fixing 
the Coulomb gauge in large volumes. We study finite-size and discretization effects. While finite- 
size effects can be controlled by the cone cut, and the ghost propagator and the Coulomb potential 
become scaling functions with the cylinder cut, the equal-time gluon propagator does not show 
scaling in the considered range of the inverse coupling constant. The ghost propagator is infrared 
enhanced. The Coulomb potential is now extended to considerably lower momenta and shows a 
more complicated infrared regime. The Coulomb string tension satisfies Zwanziger's inequality, 
but its estimate can be considered only preliminary because of the systematic Gribov effect that 
is particularly strong for the Coulomb potential. 
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1. Motivation and Introduction 



According to the Gribov-Zwanziger scenario, confinement is related to the behaviour of certain 
gauge-variant two-point functions at small momenta. This is not only true for the Landau gauge. 
For the Coulomb gauge, similar predictions [[I]] exist: the transversal component D tr of the equal- 
time gluon propagator 1 

Dff{q) = (A? (k) A%-k)) = 8 ab Uj - M[\ D tr (q 2 ) (1.1) 

should vanish in the limit q — > 0, whereas the ghost propagator G, defined by 

G" h (q) 4l (< A ™ [M-%) = S ab Gtf) , (1.2) 



is infrared enhanced 2 . The physical momentum is given by q t = (2/a) • sax{nki/L) with Fourier 
momenta kj = —L + 1 , ... ,L. M denotes the Faddeev-Popov operator whose lattice version is given 
by 



M ab 



f^eTr [{^r fc }(^ + ^)^ . (1.3, 



xy 

Moreover, Zwanziger showed that the instantaneous Coulomb potential V cou i (r = \x — y\), which 
appears in the Coulomb gauge Hamiltonian through the elimination of longitudinal gluons and 
which is related to the Faddeev-Popov operator by 

S ab V coul (q 2 ) = 1 £ (e^-V [M- 1 (-A)M- 1 ]!':) , (1.4) 

x,y 

has to linearly rise with distance in real space. This is a necessary condition to satisfy confine- 
ment usually defined through an area law for Euclidean Wilson loops. The latter translates to a 
linearly rising bound state potential VwiisonfV)- More precisely, Zwanziger's inequality ^ says that 
Vwiison (r) < — | Vcoui (r) for large distances r. Zwanziger and Cucchieri pointed out that the time- 
time component of the equal-time gluon propagator D% b (q) = 8 ab D^{q 2 ) is also related to the 
instantaneous Coulomb potential, i.e. V cou \{r) ~ (A^(x)A^(y)) might serve as an useful estimate 
of the Coulomb potential [||]. We have tested these statements by means of lattice calculations 
in SU (3) Yang-Mills theory simulated with Wilson's plaquette action. During this conference we 
learned about a similar investigation carried out by Nakagawa et al. [ffl]. 



2. Details of the simulation 

A summary of the lattices and parameters under investigation is given in Tab. [I]. We used 
two different methods to fix the gauge. The overrelaxation algorithm with CO = 1 .70 on small and 
CO = 1.60 on large lattices is compared with the simulated annealing method (with 7}^ = 0.45 and 

'Here and in the following, the norm of the spatial vector q = (gi,<72><?3) is abbreviated by q = \q\. 
2 We use the shorthand notation k-x = 2nkiXi/Li. 
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Tfinal = 0-01 as initial and final temperatures, 1500 compound sweeps at linearly decreasing tem- 
perature, each consisting of one heatbath and three microcanonical sweeps, followed by obligatory 
finalizing overrelaxation). We investigated the Gribov effect by means of the first copy (fc) - best 
copy (be) strategy that was used in Landau gauge studies Since this approach is very demand- 
ing in terms of computing time, we have restricted the study of the Gribov problem to smaller 
lattices with L 4 < 24 4 . On larger lattices only one gauge copy per configuration was generated by 
means of the simulated annealing algorithm. The bare ghost and equal-time gluon propagator are 
multiplicatively renormalized at a scale jj, 2 = 9 GeV 2 . In the case of the Coulomb potential, we use 



an ultraviolet fit to Eq. 4. 1 to fix the physical scale in y-direction. 



[GeV] a [fm] V[fm 4 ] #conf 



12 4 
16 4 
24 4 
32 4 
48 4 



12 4 



16 4 



5.8 
5.8 
5.8 
5.8 
5.8 



6.0 
6.0 
24 4 6.0 
32 4 6.0 
48 4 6.0 



12 4 6.2 

16 4 6.2 

24 4 6.2 

32 4 6.2 



1.446 
1.446 
1.446 
1.446 
1.446 



N SA 



0.1364 1.64 

0.1364 2.18 4 

0.1364 3.27 4 

0.1364 4.36 4 

0.1364 6.55 4 



100 

40 

30 

30 

20 



2.118 
2.118 
2.118 
2.118 
2.118 



1.12 



T 



1.49 4 



0.0932 
0.0932 
0.0932 2.24 4 
0.0932 2.98 4 
0.0932 4.48 4 



100 

60 

40 

30 

20 



2.914 
2.914 
2.914 
2.914 



0.0677 0.81 4 

0.0677 1.08 4 

0.0677 1.62 4 

0.0677 2.17 4 



100 
40 
30 
20 



10 
15 
20 

1 

1 



20 
30 
40 

1 

1 



10 
15 
20 
1 



N OR 



20 
30 
40 



Table 1: Ensembles and corresponding parameters as used in our investigation. and Nz? indicate the 
number of gauge copies per configuration generated with simulated annealing and overrelaxation, respec- 
tively. 



3. Gauge fixing 

Gauge fixing is performed by adjusting the gauge transformations g(x,t) G SU (3) separately 
on each time-slice t such that 

1 3 

F u is] = T7J E L Tr ( 1 - -» minimum . (3.1) 

iL x i=\ 

The measurement of the D44 propagator requires also to fix the remnant global gauge freedom with 
the temporal links involved, i. e. U(if),A ~^ St U(x,t)A 8t+i- A possible Gribov effect on D44 related 
to an ambiguity in fixing the remnant gauge freedom was not studied. Note also that minimizing 
Fu[g] on the largest lattices was faced with severe critical slowing-down. About 5% of the time- 
slices did not converge within 40000 (finalizing) overrelaxation steps. For these time-slices, the 
entire gauge fixing procedure was repeated until gauge fixing could be successfully accomplished. 
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Figure 1: The Gribov ambiguity for different lattice sizes with OR and SA gauge fixing. Upper panels: 
The convergence of the relative difference SF = 1 — F cbc /F bc between the currently best functional value 
pcbc a f ter copies and the overall best value F bc . Lower panels: the bare ghost propagator at the smallest 
lattice momentum measured for the currently best gauge copy after N cp copies. 



4. Results 

4.1 Gribov effect and its impact on the observables 

Following the fc-bc strategy, we have first checked the number of gauge copies necessary to 
find a functional value F close to the global maximum. This was done separately for both gauge 
fixing methods. For details of this strategy, we refer to analogous investigations in the Landau 
gauge The result is shown in Fig. [j] and points out that simulated annealing clearly needs less 
gauge copies to assure the convergence of the functional and hence of the observables. Therefore, 
simulated annealing should be the method of choice, in particular when going to larger lattices. 
Fig. H| illustrates the impact of the Gribov effect by comparing the observables measured on the 
first overrelaxation gauge copy and on the best out of N^p simulated annealing gauge copies. As 
already seen in Landau gauge, the infrared ghost propagator is systematically overestimated by 
about 5%. The Coulomb potential shows a similar behaviour, but now the first overrelaxation 
copy overestimates the Coulomb potential by up to 100%. The equal-time gluon propagator is less 
affected by the Gribov effect. Its transversal component shows a very weak (1%) dependence on 
the Gribov copy in the infrared region. In contrast, the time-time component seems to be under- 
estimated systematically for large momenta. 

4.2 Finite size and discretization effects 

We worked with a large number of lattice sizes and values of the inverse coupling j8 in order 
to explore finite-size and discretization effects. We found that finite-size effects are visible but can 
be efficiently removed by the cone cut (see ^ and references therein). Fig. ||] illustrates this for 
the transversal equal-time gluon propagator. Similarly, the cylinder cut minimizes discretization 
effects in the cases of the ghost propagator (see Fig. |j) and the Coulomb potential. More problem- 
atic is that we find the transversal and the time-time gluon propagator showing strong, systematic 
discretization effects. This has also been found by Nakagawa et al. [Q]. 
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12 12 

10 10 10 10 10 10 

q 2 [GeV 2 ] q 2 [GeV 2 ] 

Figure 2: Gribov effect on the propagators and the Coulomb potential: ratios of measured quantites as 
measured on the first OR copy and on the best of iV~; SA copies. 



P = 5.8, no cuts p = 5.8, cone cut 




q 2 [GeV 2 ] q 2 [GeV 2 ] 

Figure 3: The finite-size effects (left) for the transversal gluon propagator, cured (right) by the cone cut. 



4.3 Infrared behaviour of the lattice observables 

All statements in this section rely on measurements on large lattices using a single ("first") 
simulated annealing gauge copy and an appropriate combination of the cylinder and the cone cut. 
Since the equal-time gluon propagator does not show scaling for j8 = 5.8 6.2 we cannot give 
reliable statements concerning its behaviour in the infrared regime. However, the turnover of the 
transversal component for the smallest j3 = 5.8 and q 2 < 0.2 GeV 2 is in agreement with a new in- 
frared regime of the Coulomb potential observed for momenta below this value. The ghost propa- 
gator is infrared enhanced. An infrared fit of the ghost propagator to q 2 ■ G(q 2 ) ~ q~ 2a for momenta 
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Figure 4: Discretization effects (left) of the ghost dressing function, cured (right) by the cylinder cut. 
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Figure 5: Left: the ghost propagator diverging stronger than 1/q 2 in the infrared, with a fit giving a too 
small infrared exponent. Right: the Fourier-transformed Coulomb potential times q 4 , revealing the onset of 
a new infrared regime for q 2 < 0.2 GeV 2 , with a fit giving ct cou i. 



q 2 < 0.17 GeV 2 gives a = 0.19(2). This is by far too small compared to the predictions of DSE 
studies and, together with the high value of # 2 /d.o.f. ~ 12, indicates that the momentum region 
accessible in our study is not yet in the scope of the envisaged infrared power laws. We should 
also admit that our fit is clearly limited by the small number of points available in this momentum 
region. 

The Coulomb potential is infrared divergent. Plotting q" ■ V C oui(<? 2 ) reveals that a new in- 
frared regime is opened for momenta smaller than q 2 ~ 0.2 GeV 2 . After an intermediate plateau is 
reached, q 4 • V cou i(^ 2 ) decreases further with decreasing momenta. We used this momentum region 
to estimate the Coulomb string tension to a cou i = (5.0± 1.3)<7 w ii son , even higher than Zwanziger's 
inequality requires. To do so, first an ultraviolet fit 

9n 1927T 2 q 2 



■ V, 



coul \ 



for q 2 > 20 GeV , A coul = 0.9(3) GeV (4. 1) 



121 h% 2 /A c 2 oul ) 

and an infrared fit representing a sum of a constant force and a l/q 2 Coulomb law, 



<7 4 -V c 



coul I 



8 K a cou i + 4 n C q 



for q 2 < 0.16 GeV 2 , C 



17.8±5.9 



(4.2) 
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had to be performed. Note that our value of the Coulomb string tension does not take into account 
the strong Gribov effect on the Coulomb potential in the infrared. It is strictly based on data for a 
single simulated annealing copy. Hence, the Coulomb string tension might be overestimated by a 
factor of two. This would be in better agreement with values obtained by measurements of partial 
Polyakov line correlators n. We also stress that our result is in contrast to SU (2) studies that had 
suggested a cou i — CT w ii son [@]- In view of the different results for the Coulomb potential and the 
time-time component of the gluon propagator, we conclude that the latter should not be used to 
estimate the Coulomb potential. The issue of non-scaling of the equal-time gluon propagator is 
currently investigated and might be solved by a more advanced renormalization procedure. 
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